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The Takagi function t{x) is a continuous non-differentiable function introduced by Teiji Takagi 
in 1903. It has appeared in a surprising number of different mathematical contexts, including 
mathematical analysis, probability theory and number theory. This paper surveys properties 
of this function. 
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ABSTRACT 



1. Introduction 



B 

The Takagi function t(x) was introduced by T. Takagi [81] in 1903 as an example of an ev- 
erywhere non-differentiable function on [0, 1]. It can be defined on the unit interval x G [0, 1] 
in by 

O ^ 1 



.1.1 

n=0 



where {{x)) is the distance from x to the nearest integer. Takagi defined it using binary 
expansions, and showed that his definition was consistent for numbers having two binary 
expansions (dyadic rationals). The function is pictured in Figure [!} 

An immediate generalization of the Takagi function is to set, for integer r > 2, 

oo 

^rix):=Y.^ii^^^))- (1-2) 

n=0 

In 1930 van der Waerden [85] studied the function tiq and proved its non-differentiability. In 
1933 Hildebrandt |40j simplified his construction to rediscover the Takagi function. Another 
rediscovery of the Takagi function was made by of de Rham [73J in 1957. 

The Takagi function has appeared in a surprising number of different contexts, including 
analysis, probability theory and number theory. It is a prescient example of a self-similar con- 
struction. An important feature of this function is that it satisfies many self-similar functional 
equations, some of which are related to the dilation equations that appear in wavelet theory. It 
serves as a kind of exactly solvable "toy model" representing a solution to a discrete version of 
the Laplacian operator. It serves as an interesting test case for determining various measures of 
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Figure 1: Graph of the Takagi function t(x) 



irregularity of behavior of a function. It also turns out to be related to a number of interesting 
singular functions. 

This object of this paper is to survey properties of Takagi function across all these fields. 
In Section [2] we briefly review some history of work on the Takagi function. In Sections [3] 
and |4] we review basic analytic properties of the Takagi function. In Section [5] we discuss its 
connection to dynamical systems. In Section [6] we treat its Fourier transform, and in Section 
[7] its relation with Bernoulli convolutions in probability theory. Section [8] summarizes analytic 
results quantifying the local oscillatory behavior of the Takagi function at different scales, 
which imply its non-differentiabilty. Results related to number theory appear in Sections [9] 
and |10[ connecting it with binary digit sums and with the Riemann hypothesis, respectively. 



In Section 11 we describe properties of the graph of the Takagi function. In Sections |12|- 114| we 
present results concerning the level sets of the Takagi function. These include results recently 
obtained jointly with Z. Maddock ([56], [SZ])- Many of the results in these sections are based 
on properties specific to the Takagi function. This again exhibits its value as a "toy model", 
where exact calculations are possible. 

In this paper most results are given without proof. We have included proofs of several 
results not conveniently available, and some of these may not have been noted before, e.g. 



Theorem |6 . 2| and Theorem 11.4 We also draw the reader's attention to another recent survey 
of results on the Takagi function given by Allaart and Kawamura [H], which is somewhat 
complementary to this one. 



2. History 

The Takagi function was introduced during a period when the general structure of non- 
differentiable functions was being actively explored. This started from Weierstrass's discovery 
of an everywhere non-differentiable function, which he lectured on as early as 1861, but which 
was first published (and attributed to Weierstrass) by Du Bois-Reymond [T^ in 1875. Weier- 
strass's example has no finite or infinite derivative at any point. The example of Takagi is 
much simpler, and has no point with a finite derivative, but it does have some points with a 



well-defined infinite derivative, see Theorem 8.6 below. Pinkus [70] gives additional history 
on this problem. 
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In 1918 Knopp [50] studied many variants of non-differentiable functions and reviewed 
earlier work, including that of Faber [2H], and Hardy [55] . among others. He considered 
functions of the general form 

oo 

F(x):=5]a-0(6"x), (2.1) 

n=0 

where (p{x) is a given periodic continuous function of period one, for real numbers < a < 1 and 
b > 1. This general form includes the Takagi function as well as functions in the Weierstrass 
non-differentiable function family 

oo 

Wa,b{x) := ^a"cos(7r6"x). 

n=0 

Weierstrass showed this function has no finite or infinite derivative when < a < 1 and b is an 
odd integer with ab > 1 + In 1916 Hardy |36l Theorem 1.31] proved that the Weierstrass 
function has no finite derivative for real o, b satisfying < a < 1 and ab > 1. For the special 
case where = ((x)), relevant to the Takagi function, Knopp proved ([501 P- 18]) the non- 
differentiability for real < a < 1 and b being an even integer with ab > 4. The Takagi 
function has a6 = 1, so is not covered by Knopp's result; however a later result of Behrend [T3t 
Theorem III] in 1948, applies to (/>(x) = {{x)) and establishes no finite derivative for integer b 
and < a < 1, having ab > 1, with some specific exceptions. 

In the 1930's there were significant developments in probability theory, including its for- 
malization in terms of measure theory. The expression of Lebesgue measure [0, 1] in terms 
of the induced measure on the binary coefficients, reveals that these measures are given by 
independent coin fiips (Bernoulli trials). This measure density can be expressed as an infinite 
product (Bernoulli convolution), see Kac [45j for a nice treatment. In 1934, Lominicki and 
Ulam [59] studied related measures where biased coin flips are allowed. These measures are 
generally singular with respect to Lebesgue measure. In 1984 Hata and Yamaguti [38] noted 



a relation of the Takagi function to this family of measures, stated below in Theorem 7.2 

There has been much further work studying properties of Bernoulli convolutions, as well 
as more general infinite convolutions, and their associated measures. Additional motivation 
comes from work of Jessen and Wintner [44] concerning the Riemann zeta function, which is 
described at length in 1938 lecture notes of Wintner [86]. Peres, Schlag and Solomyak [69] give 
a recent progress report on Bernoulli convolutions, and Hilberdink [39J surveys connections of 
Bernoulli convolutions with analytic number theory. 

In the 1950's Georges de Rham ([7T], [72], [73]) considered self-similar constructions of 
geometric objects, again constructing a function equivalent to the Takagi function. In this 
context Kahane [46] noted an important property of the Takagi function. Similar constructions 
appear in the theory of splines, of functions iteratively constructed using control points. 

The Takagi function appeared in number theory in connection with the summatory func- 
tions of various arithmetic functions associated to binary digits. The analysis of such sums 
began with Mirsky |66j in 1949, but the connection with the Takagi function was first observed 
by Trollope |83j in 1968. It was further explained in a very infiuential paper of Delange |26j in 
1975. A reformulation of the theory in terms of Mellin transforms was given by Flajolet et al 
[ST] . These results are discussed in Section [sj 

The Takagi function can also be viewed in terms of a dynamical system, involving iterations 
of the tent map. This viewpoint was taken in 1984 by Hata and Yamaguti [3H]- Here the Takagi 
function can be seen as a kind of fractal. For further information on the fractal interpretation 
see Yamaguti, Hata and Kigami |90) . 
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In the 1990's the construction of compactly-supported wavelets led to the study of dilation 
equations, which are functional equations which linearly relate functions at two (or more) 
different scales, see Daubechies [22]. Basic results on the solution of such equations appear 
in Daubechies and the author [23], [24]. These functions can be described in terms of infinite 
products, which are generalizations of Bernoulli convolutions (see The Takagi function 

appears in this general context because it satisfies a non-homogenous dilation equation, driven 



by an auxiliary function, which is stated in Theorem 4.1 Related connections with de Rham's 
functions were observed by Berg and Kriippel [12] in 2000 and by Kriippel [55] in 2009. 

The Takagi function has appeared in additional contexts. In 1995 Frankl, Matsumoto, 
Rusza and Tokushige [32J gave a combinatorial application. For a family of /c-element sets 
of the A^-element set [N] := {1,2,..., N} and for any i < k the shadow of J- on ^-element 

sets is the set of all ^-element sets that are contained in some set in T. The Kruskal-Katona 
theorem asserts that the minimal size of ^-shadows of size m families J- of /c-element sets 
is attained by choosing the sets of J- as the first m elements in the /c-element subsets of [N] 
ordered in the co-lexicographic order. The Kruskal-Katona function gives this number, namely 

Kf{m) := m.\n{^/S.^{F) : T consists of k-element sets, = m}. 

This number is independent of the value of N , requiring only that N be sufficiently large that 
such families exist, i.e. that (^) > m. The shadow function Sk{x) is a normalized version of 
the Kruskal-Katona function, taking i = k — 1, which is given by 

k /Ik — 1 \ 

Sk{x) := -^Kt,{l(^ ^ JxJ), 0<,x<l. 

Theorem 4 of [32] states that as A; — t- oo the shadow functions Sk{x) uniformly converge to 
the Takagi function r(x). The Takagi function also appears in models of diffusion-reaction 
processes ([35J) and in basins of attraction of dynamical systems ([87j). 

3. Basic Properties: Binary Expansions 

Takagi's definition of his function was in terms of binary expansions, which we write 

x = ^^ =0.bib2b3---, each 6j G {0,1}. (3.1) 

i=i 

The binary expansion of x is unique except for dyadic rationals x = ^ , which have two possible 
expansions. For < x < 1 the distance to the nearest integer function {{x)) is 

_ ( X if < X < ^, i.e. bi = 
^^''^ -{ 1-x if l<x<l, Le.bi = l. ^"^-^^ 



For n > 0, we have 

((2"x)) 

where we use the bar-notation 
to mean complementing a bit. 



0.bn+lbn+2bn+3--- if = 

(3.3) 

^ 0.bn+lbn+2bn+3--- if &n+l = 1, 



6 = 1-6, for 6 = or 1, (3.4) 
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Definition 3.1. Let x £ [0,1] have binary expansion x = YlJLi |j = 0.5i62&3---) with each 
bj £ {0, 1}. For each j > 1 we define the following integer- valued functions. 

(1) The digit sum function Nj{x) is 

N]{x):=hi+b2 + --- + hj. (3.5) 

We also let N^{x) = j — Nj{x) count the number of O's in the first j binary digits of x. 

(2) The deficient digit function Dj{x) is given by 

Dj{x) := N^{x) - N}{x) = j - 2N]{x) = j - 2(6i + ^2 + • • • + bj). (3.6) 

Here we use the convention that x denotes a binary expansion; dyadic rationals have two 
different binary expansions, and all functions Nj{x), Nj{x), Dj{x) depend on which binary 
expansion is used. (The name "deficient digit function" reflects the fact that Dj{x) counts the 
excess of binary digits 6^ = over those with 6^ = 1 in the first j digits, i.e. it is positive if 
there are more O's than I's. ) 

Takagi's original characterization of his function, which he used to establish non-differentiability, 
is as follows. 

Theorem 3.1. (Takagi 1903) For x = 0.bib2bs--- the Takagi function is given by 

m=l 

in which < £m = £m{x) < m — 1 is the integer 

£m{x) = #{i: l<i<m, bi^bn,}. (3.8) 
In terms of the digit sum function N}^{x) = 6i + 62 + ••• + 

r Niix) if 5^+1 = 0, 

£„,+i{x) = I (3.9) 
[ m - iV^(x) if bm+i = 1- 

Dyadic rationals x = 1^ have two binary expansions, and in consequence the formulas 
above give two expansions for t{x). Theorem 3.1 asserts that these expansions give the same 
value; one may verify that t{x) itself will then be another dyadic rational, with the same or 
smaller denominator. See [SG, Lemma 2.1] for a proof of this result. 

We deduce some basic properties of the Takagi function from its definition. 

Theorem 3.2. (1) The Takagi function r(x) maps rational numbers x to rational numbers 
t{x). 

(2) The values of the Takagi function satisfy < r(x) < |. The minimal value y = is 
attained only at x = 0, 1. The maximal value 7/ = | is also attained at some rational x, in 
particular t{^) = |. 
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Proof. (1) This follows from (1.1) because for rational x the sequence ((2"a;)) takes rational 
values with bounded denominators, and becomes eventually periodic. Summing geometric 
series then gives the rationality. 

(2) The lower bound case is clear, and by inspection is attained for x = 0, 1 only. The 
upper bound is proved by checking that {{x)) + \ {{2x)) < \ holds for all x G [0, 1], and using 
this on successive pairs of terms in (1.1) to get 7"(2;)<| + | + ^ + -- - = |. One checks from 



(1.1) that r(i) = □ 



Remark. The set of values x having r(x) = | is quite large; see Theorem 12.2 Concerning 



the converse direction to Theorem 



3.2 



(1): It is not known which rationals y with < y < 



have the property that there is some rational x such that t{x) = y. 



The Takagi function t(x) can be constructed as a limit of piecewise linear approximations. 
The partial Takagi function of level n is given by: 



i=o 

See Figure[2]for T2{x),t^{x) and T4(a 



[X). 



Theorem 3.3. The piecewise linear function Tn{x) = X]j=o ^^^T^ linear on each dyadic 



interval 



(1) On each such interval Tn{x) has integer slope between —n and n given by the deficient 
digit function 

Dn{x) = N^{x) - Afi(x) = n - 2(6i + 62 + • • • + 6„), 



Here x = 0.6i62^3--- may be any interior point on the dyadic interval, and can also be an 
endpoint provided the dyadic expansion ending in 's is taken at the left endpoint ^ and that 
ending in I's is taken at the right endpoint ^tr". 

(2) The functions Tn{x) approximate the Takagi function monotonically from below 



Tl{x)<T2{x)<T^{x)<... (3.11) 

The values {rn(x) : n > 1} converge uniformly to r(x), with 

|r„(x)-T(x)|<^-i^. (3.12) 

(3) For a dyadic rational x = perfect approximation occurs at the n-th step, and 

t{x) = Tm{x), for all m > n. (3.13) 

Proof. All statements follow easily from the observation that each function fn{x) := 2" 
a piecewise linear sawtooth function, linear on dyadic intervals [2^, ^^ft], with slope having 
value +1 if the binary expansion of x has bn+i = and slope having value —1 if 6^+1 = 1- The 
inequality in (3.12) also uses the fact that max^g[o 1] r(x) = |. □ 
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Figure 2: Approximants to Takagi function: (left to right) r2(x), T3(x), T4(x). Slopes of linear 
segments are labelled on graphs. A vertex (x, y) is marked by a solid point if and only if 
X G 0.^, as defined in (13.1). 



The Takagi function itself can be directly expressed in terms of the deficient digit function. 



The relation (3.6) compared with the definition (3.9) of £m{x) yields 

m 1 



(x) 



'lf"^+^D^{x). 



Substituting this in Takagi's formula (3.7) and simplifying (noting that ^i(x) = Dq[x) = 0) 
yields the formula 



1 1 



r X 



n Dm{x) 



(3.14) 



\m=l 



4. Basic Properties: Functional Equations and Self-AfRne Rescalings 



We next recall two basic functional equations that the Takagi function satisfies. These have 
been repeatedly found, see Kairies, Darslow and Frank [JH] and Kairies |47] . 

Theorem 4.1. (1) The Takagi function satisfies two functional equations, each valid for 
< X < 1, the reflection equation 

r(x)=r(l-x), (4.1) 



and the dyadic self- similarity equation 



r( — ) = -X H — t(x). 
^2' 2 2' 



(4.2) 



(2) The Takagi function on [0,1] is the unique continuous function on [0,1] that satisfies 
both these functional equations. 



y, whence ( 1.1 ) 



Proof. (1) Here (4.1) follows directly from (1.1), since {{kx)) = {{k{l — x))) for G Z. To 
obtain (4.2), let x = 0.5i62&3"- and set y := | = O.O616263.... Then 
gives 

((2"2/))\ 



My) 



+ > - = x + t(x). 



\n=l / m=0 

(2) The uniqueness result follows by showing that the functional equations determine the 
value at all dyadic rationals. Indeed, the dyadic self-similarity equation first gives t(0) = 0, 
whence r(l) = by the reflection equation. Then t(^) = ^ by the self-similarity equation, and 
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now we iterate to get all other dyadic rationals. Since the dyadic rationals are dense, there is 
at most one continuous interpolation of the function. The fact that a continuous interpolation 
exists follows from Theorem 3.3, (This result was noted by Knuth |5H Exercise 82, solution 
p. 103]. □ 

The functional equations yield a self-afhne property of the Takagi function associated to 
shifts by dyadic rationals x = 

Theorem 4.2. For an arbitrary dyadic rational xq = then for x G ^fe^] given by 
X = xq + , there holds 

^(^0 + ^) = r{xo) + ^(r(y) + D„(a;o) y\ < y < 1, (4.3) 

That is, the graph of t{x) on is a miniature version of the tilted Takagi function 

t{x) + Dn{xQ)x, shrunk by a factor 4r and vertically shifted by t{xq). 



Proof. By Theorem 3.3 '1), we have t„(xo + ^) = t„(xo) + Dnixo) ■ Therefore, by (1.1) 
it follows that 

t{x) = rn{x) + Y,^^^^ 

j=n 

= Tn[Xo) + Dn{xo) ■ — + 

j=n 

= r{xo) + ^{Tiy) + Dn{xo)y). □ 



Theorem 



4.2 



simplifies in the special case of = ^ having D„(xo) = 0, which we call 



a balanced dyadic rational; such dyadic rationals can only occur when n = 2m is even. The 
formula (4.3) then becomes 

r(xo + |^)=r(xo) + ^, < y < 1, (4.4) 

which a shrinking of the Takagi function together with a vertical shift. Balanced dyadic 
rationals play a special role in the analysis of the Takagi function. 

As a final topic in this section, a convex function is characterized as one that satisfies the 
condition 

f( ^ + y ) < /(^) + fjy) 

The Takagi function is certainly very far from being a convex function. However one can 
establish the following approximate mid-convexity property of the Takagi function, due to 
Boros US] 



Theorem 4.3. (Boros 2008) The Takagi function is (0, ^)-midconvex. That is, it satisfies 
the bound 

x + y T{x) + T{y) 1 
r( — )< ^ + ^\x-y\. 

This result establishes the sharpness of general bounds for approximately midconvex func- 
tions established by Hazy and Pales [37] . 
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5. The Takagi function and Iteration of the Tent Map 



An alternate interpretation of the Takagi function involves iterations of the symmetric tent 
map T : [0, 1] — )• [0, 1], given by 

^(^) = 1 2-2x if i<x<l. (^-^^ 

This function extends to the whole real line by making it periodic with period 1, and is then 
given by T(x) = 2((x)), and it satisfies the functional equation T{x) = T{l — x). A key property 
concerns its behavior under iteration T°"'{x) := T{T°^"'~^\x)) (n-fold composition of functions 
is here denoted on). 

it satisfies the functional equation under composition 

r°"(x) = r(2"2;). (5.2) 

The tent map T(x) under iteration is an extremely special map. It defines under iteration a 
completely chaotic dynamical system, defined on the unit interval, whose symbolic dynamics is 
the full shift on two letters. This map has Lebesgue measure on [0, 1] as an invariant measure, 
and this measure is the maximal entropy measure among all invariant Borel measures for T{x). 



As an immediate consequence of (5.2) we have the following formula for t(x), noted by 
Hata and Yamaguti ^ in 1984. 

Theorem 5.1. (Hata and Yamaguti 1984) The Takagi function is given by 

oo ^ 
n=l 

where T°'^(x) denotes the n-th iterate of the tent map T{x). 



This result differs conceptually from the original definition (1.1) which produces powers of 
2 by a rescaling the variable in a fixed map, in that the powers of 2 in (5.3) are produced by 
iteration of a map. 

Hata and Yamaguti [38] defined generalizations of the Takagi function based on iteration 
of maps. They defined the Takagi class Ej- to consist of those functions given by 



/(x) = X;a„r°"(x). 



n=l 



where Yl'^=i \^n\ < oo viewed as members of the Banach space of continuous functions C°([0, 1]), 
under the sup norm. This class contains some piecewise smooth functions, for example 



oo ^ 

5:-T-(x)=x(l-x). 



4 

71=1 



Hata and Yamaguti showed that continuous functions in can be characterized in terms 
of their Faber-Schauder expansions in C'^([0, 1]), and that the Takagi class Et is a closed 
subspace of C([0, 1]). The standard Faber-Schauder functions, defined by Faber |HD] in 1910 
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and generalized by Schauder fTSJ pp. 48-49] in 1927, are l,x and {Si2j{x) : i > 0,0 < i < 
2^ — 1} given by dyadically shrunken and shifted tent maps 

1,2-1 \ J u 2-? 2-? 2-J+i 

These functions form a Schauder basis of the Banach space C'^([0, 1]), in the sup norm, taking 
them in the order l,x, followed by the other functions in the order: 5*4 2^ precedes S^, 2j> if 
j < j' or j = j' and i < i' . (For a discussion of Schauder bases see Megginson |62, Sec. 4.1].) 
A function f{x) in C'^([0, 1]) has Faber-Schauder expansion 

oo 2^-1 

f{x) = ao + aix + ^ ^ aijSi^2^{x), (5.4) 

with coefficients ao = /(O), a(l) = /(I) — /(O), and 

«U := /(^) - I {fi^) + /(^)) • 

Functions in the Takagi class satisfy the restriction that ao = ai = and their Faber-Schauder 
coefficients {aij} depend only on the level j. 

Hata and Yamaguti \38\ Theorem 3.3] proved the following converse result. 

Theorem 5.2. (Hata and Yamaguti 1984) A function f{x) G C°([0, 1]) belongs to the Takagi 

class Erp if and only if its Faber-Schauder expansion f{x) = aQ + aix + Y^^=Q Xlilo^ flij'S'j,2J (^) 
satisfies 

1. The coefficients oq = ai = and 

qqj = Oij for all j >0, < i < 2^~^. 

2. If we set Cj = aoj, then 

oo 

^\Cj\ < OO. 

j=0 

Note that if the conditions 1, 2 above hold, then f{x) is given by the expansion 

oo 

f{x) = Y,cjT°(^+^\x). 

j=0 

Hata and Yamaguti also viewed functions in the Takagi class as satisfying a difference 
analogue of Laplace's equation, using the scaled central second difference operators 

A,,2. (/) := /(^) + /(^) - 2/(^), 

for < i < 2-' — 1, j > 0, along with Dirichlet boundary conditions. (The Faber-Schauder 
coefficients Oij = — ^Aj2i(/)-) They obtained the following existence and uniqueness result 
([Ml Theorem 4.1]). 
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Theorem 5.3. (Hata and Yamaguti 1984) Given data {cj : j > 0}, the infinite system of 
linear equations defined for 

Ai,2.(/) = cj, J >0, 0<i<2^ -l, 

has a continuous solution fix) G C''([0, 1]) satisfying Dirichlet boundary conditions /(O) = 
/(I) = if and only if 

oo 

^\Cj\ < OO. 



j=0 



In this case f{x) £ Et, with 



^ oo 

f{x) = --Y,CjT°^^^'\x). 



2 

3=0 

There are interesting functions obtainable from the Takagi function by monotone changes of 
variable. The tent map T{x) is real-analytically conjugate on the interval [0, 1] to a particular 
logistic map 

F{y):=Ay{l-y), < y < 1. (5.5) 
That is, F{y) = ip~^ oTo (p{y) for a monotone increasing real-analytic function (/^(x), which is 

TT 

with its functional inverse 93~^(x) given by ^{x) := sin^ (^)- The dynamics under iteration of 
the logistic map (5.5 ) has been much studied. It is a post-critically finite quadratic polynomial, 
and it is affinely conjugat^to the monic centered quadratic map F{z) = — which specifies 
a boundary point of the Mandelbrot set. 

Its n-th iterate F°"'{y) is a polynomial of degree 2", which is conjugate to the Chebyshev 
polynomial r2"(x). The analytic conjugacy above gives 

^ o F°"(y) = T°" o ^{y), for n > 1. 

The change of variable x = ip{y) applied to a function f{x) = Yl'^=i CnT°'^{x) in the Takagi 
class yields the rescaled function 

oo 

9{y) = f{v{y)) = E ° (5-6) 

n=l 

motivating the further study of maps of this form. 

More generally, given a dynamical system obtained by iterating a map T : [0, 1] — )• [0, 1] of 
the interval, and a rescaling function : [0, 1] — )• C, Hata and Yamaguti define the generating 
function 

oo 

F{t,x) := ErVo^°"(^)- 

n=0 

Here these functions, where t may vary, encode various statistical information about the dis- 
crete dynamical system T. For more information on this viewpoint, see Yamaguti, Hata and 
Kigami [901 Chap. 3]. 



The conjugacy is i/ = ^(z) ~ + | with inverse z = ip ^(y) = — 4j/ + 2. 



11 



6. Fourier Series of the Takagi Function 



The Takagi function defined by (1.1) extends to a continuous periodic function on the real 
hne with period 1, which comes with a Fourier series expansion. The functional equation 
r(x) = r(l — x) for < X < 1 implies that the extended function also satisfies 

r(x) = r(— x), (6-1) 

so that it is an even function. As a consequence, its Fourier series only involves cos(27rn2;) 
terms. It is easily computable from the known Fourier series of the symmetric tent map, as 
follows. 

Theorem 6.1. The Fourier series of the Takagi function is given by 

:=^c„e2™^ (6.2) 



T(X) 



in which 



1 1 

Co = / T{x)dx = - (6.3) 







and for n > there holds 

1 



-n — '-—n 



2^(2^ + l)27r2^ 



where n = 2"(2/c + l). (6.4) 



Proof. The even function 'ip{x) = {{x)) has Fourier series with real coefficients 6„ = b-n, 
given as 



oo 



n=—oo 

1 

with b^n = J\ \x\e'^'^^'^^dx. Clearly &o = | and, integrating by parts. 



1 1 

m 

1 , 1 



2ivinx-\ \^~2 



ATTin ^(27rm)2 ^ '^=0 

\n - ^i^' " odd, 
1 



iS^' " even. 



A similar calculation on the interval [—^,0] gives the complex conjugate value, whence 

-T^-g if n is odd, 

if n 7^ is even. 



bn = 6-n = I 

Now r(x) = Ylm=o 2^((2'"^))' the Fourier coefficients of il>{2'^x) are 6™ := 6^/2™. if 2™'|n 
and otherwise. By uniform convergence of the sum we obtain cq = X^^o 2^^0 = \^ and, for 



n = ±2™(2A; + 1), we obtain 



00 ^ m ^ 
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which is the result. □ 



Note that the decay of the Fourier coefficients as n — )• oo has Hmsup„_i.oo ''^^|cn| > 0. This 
fact directly implies that t(x) cannot be a C^-function. However much more about its oscillator 
behavior, including its non-differentiability, can be proved by other methods. 

As a direct consequence of this result, the Takagi function is obtainable as the real part of 
the boundary value of a holomorphic function on the unit disk. 

Theorem 6.2. Let {c„ : n G Z} he the Fourier coefficients of the Takagi function, and define 
the power series 

^ oo 

/W = 2Co + ^c„z". (6.5) 

n=l 

This power series converges absolutely in the closed unit disk {z : \z\ < 1} to define a holo- 
morphic function in its interior, and it has the unit circle as a natural boundary to analytic 
continuation. It defines a continuous function on the boundary of the unit disk, and its values 
there satisfy 

fie'-^^) = liXi9) + ^Y{e)) (6.6) 
in which X{9) = t{9) is the Takagi function. 

Remark. The imaginary part defines a function Y[9) is a new function which we term the 
conjugate Takagi function. 

Proof. The Fourier coefficients of the Takagi function satisfy 

I 1-1 ^ 1 ^ _ 1 1 

2^ |cn| - 2 + 2^ 2— 1 lz^(2fc + l)2^2 ) - 2 + 2 

n=0 m=0 \fc=0 ^ ' / 

It follows that the power series for f{z) converges absolutely on the unit circle, so is continuous 
on the closed unit disk, and holomorphic in its interior. Since the Fourier series for the Takagi 
function is even, by inspection 

/(e^^'^) + /(e"^^^^) = 2i?e(F(e2^^^)) = X{e). 



This justifies (6.6), defining U{9) by 

/(e^^'^) - /(e~^^^^) = 2iIm{F{e^^'^)) = iY{e). 

Since the Takagi function is non-differentiable everywhere, the function f(z) cannot analyti- 
cally continue across any arc of the unit circle, so that the unit circle is a natural boundary 
for f{z). □ 

As mentioned in Section[5} the Takagi function can also be studied using its Faber-Schauder 
expansion, rather than a Fourier expansion. In 1988 Yamaguti and Kigami [91] defined for 
p > 1 the Banach spaces 



Hp■.= {fix) = Y,c^,J2-i+'S,^^,ix)}, 
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defined using a rescaled Schauder basis and the norm ||/||p = ■ They deduced 

using its Schauder expansion that the Takagi function belongs to the Banach space Hp for all 
p > 2 (It does not belong to H2 , which is a Hilbert space coinciding with a space denoted in 
j91|.) In 1989 Yamaguti |89j proposed a generalized Schauder basis consisting of polynomials 
in X, in which to study the Takagi function and other functions. 



7. The Takagi Function and Bernoulli Convolutions 

The Takagi function also appears in the analysis of non-symmetric Bernoulli convolutions. Let 

00 ^ 



i=i 



in which the digits ej £ {0, 1} are drawn as independent (non-symmetric) Bernoulli random 
variables taking value with probability a and 1 with probability 1 — a. Then a; is a random 
real number in [0, 1] with cumulative distribution function 



X 



La{x) = fia{[0,x]) = / dHa, 

Jo 

in which /Xq, is a certain Borel measure on [0, 1]. These functions were introduced in 1934 by 
Lomnicki and Ulam |59| pp. 267-269] with this interpretation. The measure fia is Lebesgue 
measure for a = ^, and is a singular measure otherwise. 

In 1943 R. Salem |74j gave a geometric construction of monotonic increasing singular func- 
tions that includes these functions special case. 

Theorem 7.1. (Salem 1943) The function La{x) has Fourier series La{x) ~ Yl'n=i 27rnx 
using the formula (for real t ) 



/ e^^''''dL^{x) = W[a + {l-a)e^y 



fc=i 

The Fourier coefficients are given by infinite products 



00 

^.•„ -i-r / vrn , , 7rn\ 

k=l 



An interesting property of the Lomnicki-Ulam function La{x) is that it satisfies a two-scale 
dilation equation 



La{x) 



aLa{2x) for < x < 

{I - a)La{2x - I) + a for ^ < rr < 1. 



In 1956 de Rham ([7T] , [72], [7H]) studied such functional equations in detail. Dilation equa- 
tions are relevant to the construction of compactly supported wavelets, and general properties 
of solutions to such equations were derived in Daubechies and Lagarias [53j, [25] . 

In 1984 Hata and Yamaguti [38, Theorem 4.6] made the following connection of the 
Lomnicki-Ulam functions La{x) to the Takagi function, cf. \38\ p. 195]. 
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Theorem 7.2. (Hata and Yamaguti (1984)) For fixed x the function {a) := Lq,(x), initially 
defined for a £ [0, 1], extends to an analytic function of a on the lens-shaped region 

-D = {a G C : |a| < 1 and \1 — a\ < 1}. 

The Takagi function appears as the derivative of these functions at the central point a = \: 

^K{x)1^^=2t{x). (7.1) 

This result permits an interpretation of the Takagi function as a generating function of 
a chaotic dynamical system (Yamaguti et al. [90< Chapter 3]). This result was generalized 
further by Sekiguchi and Shiota [76] , 



8. Oscillatory Behavior of the Takagi Function 

The main feature of the Takagi function is that it is non-differentiable everywhere. Takagi |81j 
showed the function has no two-sided finite derivatives at any point. In 1984 Cater [20] showed 
that the Takagi function has no one-sided finite derivative at any point. However it does have 
well-defined (two-sided) improper derivatives equal to +oo (resp. — oo) at some points, see 
Theorem [HH 

The non-differentiability of the Takagi function is bound up with its increasing oscillatory 
behavior as the scale decreases. Letting < h < 1 measure a scale size, we have the following 
elementary estimate bounding the maximal size of oscillations at scale h. 

Theorem 8.1. For 0<x<x + h< 1, the Takagi function satisfies 

|r(x-h/i) -r(x)| < 2/ilog2^. (8.1) 



Proof. Suppose 2 ^ < h <2 , so that n < log2 ^. Theorem 3.3 gives the estimate that 

2 

3 2" 



|r(x) — T„(x)| < 14^, and we know r„(x) has everywhere slope between —n and n. It follows 



that 

|r(x + h)- t{x)\ < ^ +nh< h{n + ^) < Ihlog^ ^. 

as required. □. 



Theorem 8.1 is sharp to within a multiplicative factor of 2, since for h = 2 

n , , 1 



r(/i)-r(0)=r(2-") = - = /ilog2^. 



In fact the multiplicative factor of 2 can be decreased to 1 as the scale decreases, in the following 
sense (Kono [52} Theorem 4]). 

Theorem 8.2. (Kono 1987) Let au{h) = log2 ^. Then there holds 

r(x)-r(y) 
iimsup 1 r- = 1. 

\x-y\^o+ \x - y\ cru{\x - y\) 

and 

t{x) - T{y) _ 



liminf 

\x-y\^o+ \x - y\ au{\x - y\) 
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The average size of extreme fluctuations for most x is of the smaUer order /i y log2 |; y 2 log log log2 ^ , 
as given in the following result. (Kono [52\ Theorem 5]). 



Theorem 8.3. (Kono 1987) Let ai{h) = Ylog2 ^. Then for almost all x G [0, 1] there holds 

T{x + h)-T{x) 

limsup = 1, 

h^o+ hai{h)^/2loglogai{h) 

and 

. p t(x + h) - t(x) 
lim mt — 



h^o+ hai{h)^/2\og\ogai{h) 
Kono used expansions of t{x) in terms of Rademacher functions to obtain his results. 



Finally, if one scales the oscillations by the factor /ly log2 ^ on the scale size h then one 
obtains a Gaussian limit distribution of the maximal oscillation sizes at scale /i as /i — )• 0^ , in 
the following sense (Gamkrelidze[3^ Theorem 1]). 



Theorem 8.4. (Gamkrelidze 1990) Let ai{h) = \ \0g2 \ - Then for each real y, 



y 



lim Meas {x : ^^"^ | /^^^ <y] = ^ I e-p^dt. (8.2) 
h^o+ hai{h) \/2tt J-00 



Gamkrelidze also observes that Theorem 8.3 can be derived as a consequence of this result. 
Namely, in this context Theorem 8.3 is interpretable as analogous to the law of the iterated 
logarithm in probability theory. 

The oscillatory behavior of the Takagi function has also been studied in the Holder sense. 
We define on [0, 1] the class of continuous functions, the class of continuously dif- 
ferentiable functions (with one-sided derivatives at the endpoints) and, for < a < 1, the 
(intermediate) Lipschitz classes 

Lzp":={/GC°: there exists > with |/(x) - /(y)| < - y|", x,yG[0,l]}. (8.3) 

Theorem 8.5. (de Vito 1985; Brown and Kozlowski 2003) 
(1) The Takagi function r belongs to the function class 

T€ f] Lip''. (8.4) 

0<a<l 

(2) The Takagi function r does not agree with any function g ^ on any set of positive 
measure. In fact, if M is a subset of [0, 1] with positive measure then the set 

D(t, M) := {iM^lM : x,yGM with x ^ y} (8.5) 
y -X 

is unbounded. 
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Proof. (1) This was shown by de Vito [83] in 1958. He showed moreover, that for < a < 1 
and < x,y < 1, there holds 



9« — 1 



|r(x) - r(y)| < 7z^\x - y 



1 



Another proof of (8.4) was given by Shidfar and Sabetfakhri \77\ . An extension to r,. for ah 



even r > 2 in (1.2) follows from results in Shidfar and Sabetfakhri |78j . 



2) A proof establishing this and (8.5) was given by Brown and Kozlowski in 2003. □ 



Concerning the non-differentiability of the Takagi function, Takagi established there is no 
finite derivative at any point of [0, 1]. However the Takagi function does have some points where 
it has a well-defined ( two-sided) infinite deriviative. These points were recently classified by 
Allaart and Kawamura [5j, who proved the following result ([5", Corollary 3.9]). 

Theorem 8.6. (Allaart and Kawamura 2010) The set of points where the Takagi function has 
a well-defined two-sided derivative r (x) = -\-oo, and the set of points where it has r (x) = — oo 
both are dense in [0, 1] and have Hausdorff dimension 1. 

There has been much further work studying oscillatory behavior of the Takagi function in 
various metrics. We mention the work of Buczolich |18j and of Allaart and Kawamura [1], 
[5]. In addition Allaart jl] has proved analogues of many results above for a wider class of 
functions in the Takagi class. 



9. The Takagi function and Binary Digit Sums 

The Takagi function appears in the analysis of binary digit sums. More generally, let Ar{x) 
denote the sum of the base r digits of all integers below x, so that 

H 

Ar{x):=Y,Sr{j). (9.1) 

j=l 

in which Sr{i) denotes the sum of the digits in the base r expansion of j. In 1949 Mirsky |66j 
observed that ^ 

Ar{x) = -{r — l)x logj. X — Er{x) 

with remainder term Er{x) = 0{x). In 1968 Trollope |83j observed that, for r = 2, this 
remainder term has an explicit formula given in terms of the Takagi function. 

Theorem 9.1. (Trollope 1968) If the integer n satisfies 2"^ < n < 2"^^^ , and if we write 
n = 2"^(1 + x) for a rational number < a; < 1, with x = ^ — 1, then 

E2{n) = 2™-^ ((1 + x) log2(l + x)-2x + r(x)) . 

This was generalized in a very influential paper of Delange f26] , who observed the following 
result. 
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Theorem 9.2. (Delange 1975) One can write the rescaled error term Ej.{n) for base r digit 
sums at integer values n as 



-Erin) = \ log^ n + Fr (log^ n) , 
n 2 

in which Fr{t) is a continuous function which is periodic of period 1. The Fourier series of the 
function Fr{t) is given by 

Fr{t) = Y,Ck{r)e''^"'\ 

The Fourier coefficients Ck{r) involve the values of the Riemann zeta function ^(s) at the points 
on the imaginary axis. 

In view of Trollope's result, in the case r = 2 the Takagi function t(x) appears in the 
function F2{t), namely this function is 

F2(t) = 1(1 + [tj - t) + 2i+W-V(2*-W-i). (9.2) 

Also for r = 2, and k ^ 0, the Fourier coefficients of F2{t) are 

1 + ^ log 2 

and for A; = 0, 

1 3 

Co = (log27r — 1) . 

21og2^ ^ ^ 4 

Results on k-th powers of binary digit sums were given in 1977 by Stolarsky [80j , who also 
provides a summary of earlier literature. A connection of these sums with a series of more 
complicated oscillatory functions was obtained in 1986 in Coquet [21]. Further work was done 
by Okada, Sekiguchi and Shiota |67| in 1995. Recently Kriippel |I54| gave another derivation 
of these results and further extended them. 

For a related treatment of similar functions, using Mellin transforms, see Flajolet et al 
|31j . A general survey of dynamical systems in numeration, including these topics, is given by 
Barat, Berthe, Liardet and Thuswaldner [11] . 



10. The Takagi Function and the Riemann Hypothesis 

The Riemann hypothesis can be formulated purely in terms of the Takagi function, as was 
observed by Kanemitsu and Yoshimoto |49l Corollary to Theorem 5]. Their encoding of the 
Riemann hypothesis concerns the values of the function restricted to the Farey fractions Fn. 

The Farey series of level consists of all reduced rational fractions < | < 1 having 
denominator at most N, their number being = + 0{NlogN). The connection 

of the approximately uniform spacing of Farey fractions and the Riemann hypothesis starts 
with Franel's theorem, cf. Franel [33J, proved in 1924. (cf. Huxley j4H p. 36]). Since that 
time many variants of Franel's result have been established. A particular version of Franel's 
theorem that is relevant here was given by Mikolas [63], |64) . [65| . 

Uniformity of distribution of a set of points (here the Farey fractions) can be measured in 
terms of the efficiency of numerical integration of functions on the unit interval obtained by 
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sampling their values at these points. This is the framework taken in the result of Kanemitsu 
and Yoshimoto. The particular interest of their result is that it applies to the numerical 
integration of continuous functions that are not necessarily differentiable, but which belong 
to the Lipschitz class Lip^~'^ for each e > 0. The role of the Takagi function here is to be a 
particularly interesting example to which their general theorem applies. 

Theorem 10.1. (Kanemitsu and Yoshimoto 2000) The Riemann hypothesis is equivalent to 
the statement that for each e > there holds, for the Farey sequence Fn, 

V t{p) - \Fn\ [ T{x)dx = O (m+') (10.1) 

as N ^ CO. 



Proof. This result was announced in Kanemitsu and Yoshimoto [49', Corollary to Theorem 
5], and a detailed proof was given in 2006 in Balasubramanian, Kanemitsu and Yoshimoto [lOj . 
According to the Principle stated in p. 4], the exponent | + e in the remainder term in 
( 10.1 ) depends on the fact that the Takagi function is in the Lipschitz class Lip'^ for any a < 1, 
see Theorem 12.4 below. (Note that in |10l Sec. 2.2] their Takagi function T[x) is defined by 
a Fourier series, which when compared with Theorem |6.1| suggests it needs a rescaling plus a 
constant term added to agree with t{x)] this does not affect the general argument.) □ 



An interesting feature of this result is that for each the left side of (10.1) is a rational 
number. This follows since the Takagi function takes rational values at rational numbers, and 
because with our scaling of the Takagi function we have T{x)dx = ^. 



11. Graph of the Takagi Function 

We next consider properties of the graph of the Takagi function 

g(r) := {ix,T{x)) : < x < 1}. 
The local extreme points of the graph of the Takagi function were determined by Kahane 

Theorem 11.1. (Kahane 1959) 

(1) The set of local minima of the Takagi function are exactly the set of all dyadic rational 
numbers in [0, 1]. 

(2) The set of local maxima of the Takagi function are exactly those points x such that the 
binary expansion of x have deficient digit function D2n{x) = for all sufficiently large n. 



Proof. This is shown in [46, Sec. 1]. Kahane states condition (2) as requiring that the binary 
expansion of x have 6i + &2 + ■ ■ ■ + &2n = n holding for all sufficiently large n. □ 

There are uncountably many x that satisfy condition (2) above. It is easy to deduce from 
this result that the graph of the Takagi function contains a dense set of local minima and local 
maxima, viewed from either of the abscissa or ordinate directions. Allaaert and Kawamura |1] 
further study the extreme values of functions related to the Takagi function. 

The Hausdorff dimension of the graph of the Takagi function was determined by Mauldin 
and Williams (61] . 

Theorem 11.2. (Mauldin and Williams 1986) The graph Q{t) of the Takagi function has 
Hausdorff dimension 1 as a subset ofM.'^. 
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Proof. This result is given as \61\ Theorem 7]. □ 

Mauldin and Williams [61] raised the question of whether the graph of the Takagi function 
is cj-finite. This was answered in the affirmative by Anderson and Pitt fS]. 

Theorem 11.3. (Anderson and Pitt 1989) The graph Q{t) of the Takagi function has a -finite 
linear measure. 



Proof. This result is given as [8", Thm. 6.4, and Remark p. 588]. □ 

In the opposite direction, one can show the graph of the Takagi function has infinite length. 
We have the following result. 

Theorem 11.4. The graph G{t) has infinite length locally. That it, it has infinite length over 
any nonempty open interval xi < x < X2 in [0, 1] . 



Proof. This is proved using the piecewise linear approximations t„(x) to the Takagi function. 
This function has control points at x = ^ where it takes values agreeing with r(x). It suffices 
to show that the length of r„(x) becomes unbounded as n — )• oo. This shows the whole graph 
has infinite length. The self-similar functional equation then implies that any little piece of 
the graph over the interval xi < x < X2 also has infinite length. 

To show the unboundedness of the length of Tn{x), as n — )• oo, it suffices to show that the 
average size of the (absolute value) of the slope of r„(x) becomes unbounded. To do this we 
note that the subdivision from level n to level n + 1 replaces each slope m with two intervals of 
slopes m + l and m — 1. These do not change the average value of the (absolute value) of slope, 
except when m = 0. Such slopes occur only on even levels 2n and there are ^^^^ of them 
(they are the balanced dyadic rationals at that level). Consequently the average value of the 
slope at level 2n (resp. 2n + 1) obeys the recursion: a2n = a2n-i and a2n+i = cL2n + 2^ 
Since 

1 /2n\ 1 



iTrn 



and y'^Li — diverges, we have a„ — t- oo as n — )• oo, proving the result. □. 



In 1997 Tricot [82] defined the notion of irregularity degree of a function / to be the 
Hausdorff dimension of its graph Q{f). He introduces a two-parameter family of norms A"'^ 
to measure oscillatory behavior. As an example, he computes these values for the Takagi 
function and its relatives ([82, Cor. 6.5, Cor. 6.7]) 



12. Level Sets of the Takagi Function 

We next consider properties of level sets of the Takagi function. For y G [0, |] we denote the 
level set at level y by 

L{y) := {x : r(x) = y, < x < 1}. 
The level sets have a complicated and interesting structure. 

Theorem 12.1. (1) The Takagi function has level sets L{y) that take a finite number of 
values, resp. a countably infinite set of values, resp. an uncountably infinite set of values. 
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Specificially, the level set L{0) = {0, 1} takes two values, the level set L{^) takes a countably 
infinite set of values, and the level set L('^) takes an uncountably infinite set of values. 

(2) The set of levels y such that L{y) is finite, resp. countably infinite, resp uncountably 
infinite, are each dense in [0, |]. 

Proof. (1) It is clear that L{0) = {0, 1} is a finite level set. A finite example given by Knuth 
[5T] is 1(1) = {x,l- x} with X = Knuth [SU pp. 20-21, 32-33] also showed that certain 

dyadic rationals x = ^ have countably infinite level sets; these include x = ^. The author 
and Maddock |56, Theorem 7.1] also show that L{^) is countably infinite. For uncountably 
infinite level sets, Baba [Qj showed that L{'^) has positive Hausdorff dimension, which implies 
it is uncountable. 

(2) For uncountably infinite sets, this follows from the self-similarity relation in Theorem 



4.2, For countably infinite case it follows from Knuth's results, see also Allaart f3]. For the 



finite case, it follows from Theorem 12.4 below. □. 



We next consider bounds on the Hausdorff dimension of level sets. 

Theorem 12.2. (Baba 1985; de Amo et al 2011) 

(1) The maximal level set L{^) has Hausdorff dimension ^. The set of levels having level 
set of Hausdorff dimension ^ is dense in [0, |]. 

(2) Every level set L{y) of the Takagi function has Hausdorff dimension at most |. 
Proof. (1) In 1984 Baba [9j showed that has Hausdorff dimension |. Using the self- 



similarity Theorem 4.2 for balanced dyadic rationals, the same holds for y = t{x) such that 
X has a binary expansion whose purely periodic part is (01)°°, and whose preperiodic part has 
an equal number of zeros and ones. The set of such x is dense in [0, 1], so their image values 
y are dense in [0, |]. 

(2) An upper bound on the Hausdorff dimension of 0.699 was established in 2010 by Mad- 
dock [60j. Recently de Amo, Bhouri, Dfaz Carrillo and Fernandex- Sanchez ^ improved on his 
argument to establish the optimal upper bound of ^ on the Hausdorff dimension. □. 

We next consider the nature of "generic level sets" , considered in two different senses. The 
first is, to draw a abscissa value x at random in [0, 1] with respect to Lebesgue measure, and 
to ask about the nature of the level set L(r(x)). The second is, to draw an ordinate value y 
at random in [0, |], and ask what is the nature of the level set L{y). Figure [s] illustrates the 
two senses. 

For sampling random abscissa values the result is as follows. 

Theorem 12.3. For a full Lebesgue measure set of abscissa points x £ [0, 1] the level set 
L(r(x)) contains an uncountable Cantor set. 

Proof. This is a corollary of a result of the author and Maddock \56\ Theorem 1.4], The 
latter result concerns local level sets, discussed below in Section fTSl □. 



For sampling random ordinate values the corresponding result is quite different. 

Theorem 12.4. (Buczolich 2008) For a full Lebesgue measure set of ordinate points y € [0, |] 
the level set L{y) is a finite set. 
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Figure 3: Ordinate level set L{y) at y = 0.5 and abscissa level set L(r(x)) at x = 0.3. 



Proof. This result was proved by Buczolich [19j in 2008. □ 

A refinement of the ordinates result is as follows. 
Theorem 12.5. Letting \L{y)\ denote the number of elements in L[y), if finite, and setting 

2 

= otherwise, one has Jq \L{y)\dy = +oo. That is, the expected number of elements in 



l-^(y)l for y G [0, i] drawn uniformly, is infinite. 



Proof. This result was obtained in |56| Theorem 6.3]. A simplified proof was given by Allaart 

m- □ 

The ordinate and abscissa results on the size of level sets are not contradictory; sampling 
a point X on the abscissa will tend to pick level sets which are "large" . In order for these two 
results to hold it appears necessary that the Takagi function must (in some imprecise sense) 
have "infinite slope" over part of its domain. 

The following result goes some way towards reconciling these two results, by showing that 
the set of large "ordinate" level sets has full Hausdorff dimension. 

Theorem 12.6. Let Tf^'^ be set of ordinates y G [0, |] such that the Takagi function level set 
L{y) has positive Hausdorff dimension, i.e. 

T'S'' := {y : dimH(L(y)) > 0}. 

Then F^'^ has full Hausdorff dimension, i.e. 

dimniT'S'') = 1. (12.1) 

It also has Lebesgue measure zero. 

Proof. The Hausdorff dimension 1 property is shown by the author and Maddock [57^ The- 
orem 1.5]). The Lebesgue measure zero property follows from Theorem 12.4[ 1). □. 



There is potentially a multifractal formalism connected to the Hausdoff dimensions of the 
sets F^'^(a) := {y : dim//(L(y)) > a}, see |57| Sec. 1.2], and for multifractal formalism see 
Jaffard [42], [13]. 

Theorem 112.41 asserts that almost all level sets in the ordinate sense are finite. The finite 
level sets have an intricate structure, which is analyzed by Allaart [3]. Finite level sets must 
have even cardinality, and all even values occur. 



22 



13. Local Level Sets of the Takagi Function 



The notion of local level set of the Takagi function was recently introduced by the author and 
Maddock [56j. These sets are closed subsets of level sets that are directly constructible from 
the binary expansions of any one of their members. 

Local level sets are defined by equivalence relation on elements x G [0, 1] based on properties 
of their binary expansion: x = • • •• Recall that in Section [2] we defined the deficient 

digit function 

Dj{x) ■.= j -2{bi + b2 + --- + bj). 

The quantity Dj{x) counts the excess of binary digits 6^ = over those with 6^ = 1 in the 
first J digits. We then associate to any x the sequence of "breakpoints" j at which tie-values 
Dj{x) = occur, setting 

Z{x):={cm: Dc^{x)=Q]. 

where we define cq = cq{x) = and set co(x) < ci(x) < C2(x) < .... This sequence of tie- values 
may be finite or infinite, and if it is finite, ending in Cn{x), we make the convention to adjoin 
a final "breakpoint" Cn-\-i{x) = +00. We call a "block" a set of digits between two consecutive 
tie- values, 

Bk{x) := {bj : Cfc(x) < j < Cfc+i(x)}. 

Two blocks are called equivalent, written Bk{x) ~ Bk'{x'), if their endpoints agree (cfc(x) = 
Cfc/(x') and Ck+i{x) = Ck'+i{x')) and either Bk{x) = By{x') or Bk{x) = Bki{x'), where the bar 
operation flips all the digits in the block, i.e. 

bj h-> 6^. := 1 - bj, Ck < j < Cfe+i. 

Finally, we define the equivalence relation on binary expansions x ^ x' to mean that Z{x) = 
Z{x'), and furthermore every block Bk{x) ~ Bk{x') for A; > 0. We define the local level set Lf'°'^ 
associated to x to be the set of equivalent points, 

4°":={x': x'^x]. 



It is easy to show using Takagi's formula (Theorem 3.1) that the relation x ~ x' implies that 
t{x) = t{x') so that X and x' are in the same level set of the Takagi function. 

Each local level set L^°'^ is a closed set, and is either a finite set if Z{x) is finite, or else is 
a perfect totally disconnected set (Cantor set) if Z{x) is infinite. The case of dyadic rationals 
X = ^ is exceptional, since they have two binary expansions, and we remove this ambiguity 
by taking the binary expansion for x that ends in zeros. 

The definition implies that each level set L{y) partitions into a disjoint union of local level 
sets A priori this union may be finite, countable or uncountable. Finite and countable 

examples are given by the author and Maddock [57] . For the uncountable case, see Theorem 



13.5 below. Note that all countably infinite level sets necessarily are a countable union of finite 
local level sets. The only level sets currently known to be countably infinite are certain dyadic 
rational levels, including x = ^. 

The paper [SB] characterizes the size of a "random" local level set sampled by randomly 



drawing an abscissa value x, as follows; this result immediately implies Theorem 12.1 



Theorem 13.1. For a full Lebesgue measure set of abscissa points x G [0, 1] the local level set 
L^°'^ is a Cantor set of Hausdorff dimension 0. 



23 



Proof. This is shown in ^56j Theorem 1.4]. □. 

In order to analyze local level sets in the ordinate space < y < |, we label each local 
level set by its leftmost endpoint. We define the deficient digit set by the condition 

:= {x G [0, 1] : x = 0.&i62&3---such that Dj{x) > 0, j = 1, 2, 3, ...} (13.1) 

This set turns out to be quite important for understanding the Takagi function. 

Theorem 13.2. (1) The deficient digit set Q,^ is the set of leftmost endpoints of all local level 
sets. 

(2) The set VL^ is a closed, perfect set (Cantor set). It has Lebesgue measure 0. 

(3) The set $7^ has Hausdorff dimension 1. 

Proof. (1) and (2) are shown in \5()\ Theorem 4.6]. 
(3) This was shown in |57l Theorem 6.1]. □ 

We show that the Takagi function restricted to the set is quite nicely behaved, as given 
in the following result. 

Theorem 13.3. The function (x) defined hyT^{x) = t{x)+x for x G 0^ is a nondecreasing 
function on il^. Define its extension to all x G [0, 1] by 

T^{x) := sup{r'^(xi) : xi < x with xi G Q^}. 

Then the function t^{x) is a monotone singular function. That is, it is a nondecreasing 
continuous function having t'^(O) = 0,t^{1) = 1, which has derivative zero at (Lebesgue) 
almost all points of [0,1]. The closure of the set of points of increase ofT^{x) is the deficient 
digit set . 

Proof. This is shown in |56| Theorem 1.5]. □. 

1,0 
0,8 
0,6 
0,4 
02 




Figure 4: Graph of Takagi singular function t^{x) 



We call the function t^{x) : [0, 1] — t- [0, 1] the Takagi singular function; it is pictured in 
Figure |4 Using the functional equation t{^x) = |(r(a;) + x) given in Theorem 4.1 '1), we 
deduce that the function ^t^{x) agrees with t(x) on the set ^0.^. This shows that the Takagi 
function is strictly increasing when restricted to the domain 

Associated to the Takagi singular function is a nonnegative Radon measure d^s, which we 
call the Takagi singular measure, such that 



r^(x) 



dfis- 



(13.2) 
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It is singular with respect to Lebesgue measure and defines an interesting probabiHty measure 
on [0,1]. The Takagi singular measure is not translation-invariant, but it has certain self- 
similarity properties under dyadic rescalings, compatible with the functional equations of the 
Takagi function. These are useful in explicitly computing the measure of various interesting 
subsets of ft^, see [57J. One may compare its properties with those of the Cantor function, see 
Dovghoshey et al [271 Sect. 5]. A major difference is that its support Q,^ has full Hausdorff 
dimension, while the Cantor set has Hausdorff dimension log3(2) ~ 0.63092. 

By an application of the co-area formula of geometric measure theory for functions of 
bounded variation (in the version in Leoni |58l Theorem 7.2 and Theorem 13.25]) to a relative 
of the Takagi singular function, the author with Maddock [5^ determined the average number 
of local level sets on a random level. 

Theorem 13.4. For a full Lebesgue measure set of ordinate points y G [0, |] the number 
N^°'^{y) of local level sets at level y is finite. Furthermore 



N'"'{y)dy = 1. 







That is, the expected number of local level sets at a uniformly drawn random level in [0, 3] is 
exactly |. 

Proof. This is shown in [56j Theorem 6.3]. A simplified derivation is given in [2J, that avoids 
the co-area formula □ 



One can easily derive both parts of Theorem 12.4 from this result. This result fails to give 
any information about the multiplicity of local level sets on those levels having an uncountable 
level set, because the set of such levels y has Lebesgue measure 0. 

Recently Allaart [2J established the following result about the multiplicities of local level 
sets in a level set. 

Theorem 13.5. (Allaart 2011) There exist levels y such that the level set L{y) contains un- 
countably many distinct local level sets. The set of such levels is dense in [0, |]. 

Allaart ([2] j |3]) has obtained further information on the structure of level sets and local 
level sets, classifying the values giving different types of level sets in terms of the Borel hierarchy 
of descriptive set theory. Some of these sets are complicated enough that they apparently live 
in the third level of the Borel hierarchy. 



14. Level Sets at Rational Levels 



In 2005 Knuth |5H Problem 83, p. 32] raised the question of determining which rational levels 
y = -g have level sets have L{'^) that are uncountable. 

The author with Maddock [56] answered the much easier question of determining when 
certain rational numbers x give uncountable local level sets L!'"'^ C L{t{x)). 

Theorem 14.1. (1) A rational number x = ^ g [0, 1] has an uncountable local level set L'p^ 

if and only if its binary expansion has a pre-periodic part with an equal number of zeros and 
ones, and if also its purely periodic part has an equal number of zeros and ones. 

(2) If a rational x has L^'"'^ uncountable, then L^f!'^ contains a countably infinite set of 
rational numbers. 
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This criterion imphes that dyadic rationals x = p- must have finite local level sets. For 
such values y = t(^) is also a dyadic rational. Concerning dyadic rational levels, Allaart [3] 
obtains the following much stronger result. 

Theorem 14.2. (Allaart 2011) Let y = he a dyadic rational with < y < |. Then the 
level set L{y) is either finite or countably infinite. 

The examples y = and 2/ = 5 given above show that both alternatives in this result occur. 
Furthermore all elements on the given dyadic rational level are rational. However if the set is 
countably infinite then these elements need not all be dyadic rationals. For example t{^) = ^, 
cf. \56\ Theorem 7.1]. 



15. Open Problems 

There remain many open questions about the Takagi function. We mention a few of them. 

1. Determine analytic and other properties of the conjugate Takagi function Y{x) defined 
in Theorem 16.21 

2. Consider the set of all x such that the abscissa level set L(r(x)) has Hausdorff dimension 
zero. Does this set have full Lebesgue measure in [0, 1]? 

3. Is there a continuous function on the interval for which the Takagi singular measure is a 
"natural" invariant probability measure? 

4. Determine the dimension spectrum function 

fr{a) := dim//({y : dimH{L{y)) > a}. 
(Some bounds on /r(a) are given in [S7].) 

5. Knuth's problem 7.2.1.3, Prob. 83]): Find a direct characterization of the set of 
rational y which have an uncountable level set. 



Acknowledgments. D. E. Knuth originally brought problems on the Takagi function to my 
attention via postcard. I thank S.T. Kuroda for helpful remarks on Takagi's original work, and 
J. -P. AUouche for some corrections. I thank Pieter Allaart and Zachary Maddock for helpful 
comments and for supplying additional references. I gratefully thank the reviewer for detailed 
comments and corrections. 



References 

[1] P. C. Allaart, On a flexible class of continuous functions with uniform local structure, J. 
Math. Soc. Japan 61 (2009), No. 1, 237-262. 

[2] P. C. Allaart, How large are the level sets of the Takagi function? eprint: 
|arXiv:1102.1616| vl, vl, 8 Feb 2011. 

[3] P. C. Allaart, The finite cardinalites of the level sets of the Takagi function, eprint: 
|arXiv:1107.0712[ vl, 4 Jul 2011. 



26 



[4] P. C. Allaart and K. Kawamura, Extreme values of some continuous nowhere differentiable 
functions, Math. Proc. Camb. Phil. Soc. 140 (2006), No. 2, 269-295. 

[5] P. C. Allaart and K. Kawamura, The improper infinite derivatives of Takagi's nowhere- 
differentiable function, J. Math. Anal. Appl. 372 (2010), 656-665. 



[6] P. C. Allaart and K. Kawamura, The Takagi function: a survey, ^arXiv: 1110. 1691 Oc- 
tober 8, 2011. 

[7] E. de Amo, I. Bhouri, M. Di'az Carrillo, and J. Fernandex- Sanchez, Hausdorff dimension 
of Takagi's function level sets, Nonlinear Analysis 74 (2011), No. 15, 5081-5087. 

[8] J. M. Anderson and L. D. Pitt, Probabilistic behavior of functions in the Zygmund spaces 
A* and A*, Proc. London Math. Soc. 59 (1989), No. 3, 558-592. 

[9] Y. Baba, On maxima of Takagi-van der Waerden functions, Proceedings of the American 
Mathematical Society, 91 (1984), No. 3, pp. 373-376. 

[10] R. Balasubramanian, S. Kanemitsu, and M. Yoshimoto, Euler products, Farey series and 
the Riemann hypothesis II, Publ. Math. Debrecen 69 (2006), no. 1-2, 1-16. 

[11] G. Barat, V. Berthe, P. Liardet, and J. Thuswaldner, Dynamical directions in numeration, 
Ann. Inst. Fourier 56 (2006), no. 7, 1987-2092. 

[12] L. Berg and M. Kriippel, de Rham's singular function, two-scale difference equations and 
Appell polynomials. Results Math. 38 (2000), 18-47. 

[13] F. A. Behrend, Some remarks on the construction of continuous non- differentiable func- 
tions, Proc. London Math. Soc. Series 2, (1948), 463-481. 

[14] P. Billingsley, Van der Waerden's continuous nowhere differentiable function, Amer. Math. 
Monthly 91 (1984), 373-376. 

[15] P. du Bois-Reymond, Versuch einer Classification der willkiirlichen Functionen reeller 
Argumente nach ihren Aenderungen in den kleinsten Intervallen, J. Reine Angew. Math. 
79 (1875), 21-37. 

[16] Z. Boros, An inequality for the Takagi function, Math. Inequalities & Applications 11 
(2008), No. 4, 757-765. 

[17] J. B. Brown and G. Kozlowski, Smooth Interpolation, Holder continuity and the Takagi- 
van der Waerden function, Amer. Math. Monthly 110 (2003), No.2, 142-147. 

[18] Z. Buczolich, Micro tangent sets of continuous functions, Math. Bohemica 128 (2003), 
No. 2, 147-167. 

[19] Z. Buczolich, Irregular 1-sets on the graphs of continuous functions. Acta Math. Hungar. 
121 (2008), 371-393. 

[20] F. S. Cater, On van der Waerden's nowhere differentiable function, Amer. Math. Monthly 
91 (1984), 307-308. 

[21] J. Coquet, Power sums of digital sums, J. Number Theory 22 (1986), 161-176. 



27 



[22] I. Daubechies, Ten Lectures on Wavelets, CBMS-NSF Regional Conference Series in Ap- 
plied Mathematics, 61. Society for Industrial and Applied Mathematics (SIAM), Philadel- 
pha, PA 1992. 

[23] I. Daubechies and J. C. Lagarias, Two-scale difference equations. I. Existence and global 
regularity of solutions, SIAM J. Math. Anal. 22 (1991), no. 5, 1388-1410. 

[24] I. Daubechies and J. C. Lagarias, Two-scale difference equations. II. Local regularity, 
infinite products of matrices and fractals, SIAM J. Math. Anal. 23 (1992), no. 4, 1031- 
1079. 

[25] I. Daubechies and J. C. Lagarias, Sets of matrices all infinite products of which con- 
verge, Lin. Alg. Appl. 161 (1992), 227-263. [Corrigendum/addendum, Lin. Alg. Appl. 
327 (2001), 69-83.] 

[26] H. Delange, Sur la fonction sommatoire de la fonction "somme des chiffres", Enseign. 
Math. 21 (1975), 31-47. 

[27] O. Dovgoshey, O. Martio, V. Ryazanov, M. Vuorinen, The Cantor function, Expo. Math. 
24 (2006), 1-37. 

[28] G. Faber, Uber stetige Funktionen, Math. Annalen. 66 (1908), 81-94. 

[29] G. Faber, Uher stetige Funktionen (Zweite Ahhandlung), Math. Annalen 69 (1910), No. 
3, 372-443. 

[30] G. Faber, Uher die Orthogonalfunktionen des Herrn Haar, Jahrsebericht Deutsch. Math. 
Verein. 19 (1910), pp. 104-112. 

[31] P. Flajolet, P. Grabner, P. Kirschenhofer, H. Prodingcr and R. F. Tichy, Mellin transforms 
and asymptotics: digital sums, Theor. Comput. Sci. 123 (1994), 191-214. 

[32] P. Prankl, M. Matsumoto, I. Z. Rusza and N. Tokushige, Minimum shadows in uniform 
hypergraphs and a generalization of the Takagi function, J. Comb. Theory A 69 (1995), 
125-148. 

[33] J. Franel, Les suites de Farey et les problemes des nombres premiers, Nachr. Ges. Wiss. 
Gottingen 1924, 198-201. 

[34] N. G. Gamkrclidzc, On a probabilistic property of Takagi 's function, J. Math. Kyoto Univ. 
30 (1990), No. 2, 227-229. 

[35] P. Gaspard and R. Klages, Chaotic and fractal properties of deterministic diffusion- 
reaction processes, Chaos 8 (1998), No. 2, 409-423. 

[36] G. H. Hardy, Weierstrass 's non-differentiable function. Trans. Amer. Math. Soc. 17 
(1916), no. 3, 301-325. 

[37] A. Hazy and Z. Pales, On approximately midconvex functions. Bull. London Math. Soc. 
36 (2004), no. 3, 339-350. 

[38] M. Hata and M. Yamaguti, The Takagi function and its generalization, Japan J. Appl. 
Math. 1 (1984), 183-199. 



28 



[39] T. Hilberdink, Some connections between Bernoulli convolutions and analytic number 
theory, in: Fractal gometry and applications: a jubilee of Benoit Mandelbrot, Part I, pp. 
233-271, Proc. Symp. Pure Math., No. 72, Part 1, AMS: Providence, RI 2004. 

[40] T. H. Hildebrandt, A simple continuous function with a finite derivative at no point, Amer. 
Math. Monthly 40 (1933), 547-548. 

[41] M. Huxley, The Distribution of Prime Numbers. Large sieves and zero-density theorems. 
Oxford University Press, Oxford 1972. 

[42] S. Jaffard, Multifractal formalism for functions part I: results valid for all functions, SIAM 
J. Math. Anal. 28 (1997), No. 4, 944-970. 

[43] S. Jaffard, Multifractal formalism for functions part II: Self-similar functions, SIAM J. 
Math. Anal. 28 (1997), No. 4, 971-998. 

[44] B. Jessen and A. Wintner, Distribution functions and the Riemann zeta-function. Trans. 
Amer. Math. Soc. 38 (1935), 48-88. 

[45] M. Kac, Statistical Independence in Probability, Analysis and Number Theory, Carus 
Mathematical Monograph No. 12. Math. Assoc. of America: New York, 1959. 

[46] J. -P. Kahane, Sur L'exemple donne par M. de Rham, d'une fonction continue sans derivee, 
Enseign. Math. 5 (1959), 53-57. 

[47] H.-H. Kairies, Takagi's function and its functional equations, Rocz. Nauk. Dydakt. Pr. 
Mat. 15 (1998), 73-83. 

[48] H.-H. Kairies, W. F. Darslow and M. J. Frank, Functional equations for a function of van 
der Waerden type, Rad. Mat. 4 (1988), no. 2, 361-374. (Errata: Rad Math. 5 (1989), no. 
1, 179-180.) 

[49] S. Kanemitsu and M. Yoshimoto, Euler products, Farey series and the Riemann hypothesis, 
Publ. Math. Debrecen 56 (2000), No. 3-4, 431-449. 

[50] K. Knopp, Fin einfaches Verfahren zur Bildung stetiger nirgends differenzierbarer Funk- 
tionen. Math. Z. 2 (1918), 1-26. 

[51] D. E. Knuth, The Art of Computer Programming, Vol. 4; Fascicle 3. Generating all 
Combinations and Partitions, Addison- Wesley: Upper Saddle River, NJ 2005. 

[52] N. Kono, On generalized Takagi functions. Acta Math. Hung. 49 (1987), No 3-4, 315-324. 

[53] M. Kriippel, On the extrema and improper derivatives of Takagi's continuous nowhere 
differentiable function, Rostock Math. KoUoq. 62 (2007), 41-59. 

[54] M. Kriippel, Takagi's continuous nowhere differentiable functions and binary digit sums, 
Rostock Math. Kolloq. 63 (2008), 37-54. 

[55] M. Kriippel, De Rham's singular function, its partial derivatives with respect to the pa- 
rameter and binary digital sums, Rostock Math. Kolloq. 64 (2009), 57-74. 

[56] J. C. Lagarias and Z. Maddock, Level sets of the Takagi function: Local level sets, eprint: 
[arXiv: 1009. 0855, v5, 29 July 2011 



29 



[57] J. C. Lagarias and Z. Maddock, Level sets of the Takagi function: Generic level sets, 
eprint: arXiv: 1011.3183 v4, 26 July 2011 

[58] G. Leoni, A First Course in Sobolev Spaces, Graduate Studies in Mathematics 105, AMS: 
Providence, RI 2009. 

[59] E. Lomnicki and S. M. Ulam, Sur la theorie de la mesure dans les espaces combinatoires 
et son application au calcul des probabilites I. Variables independantes, Fund. Math. 23 
(1934), 237-278. 

[60] Z. Maddock, Level sets of the Takagi function: Hausdorff dimension, Monatshefte fiir 
Math. 160 (2010), 167-186. 

[61] R. D. Mauldin and S. C. Wilhams, On the Hausdorff dimension of some graphs, Trans. 
Amer. Math. Soc. 298, No. 2 (1986), 793-803. 

[62] R. E. Megginson, An Introduction to Banach Space Theory, Graduate Texts in Math. No. 
183, Springer- Verlag: New York 1998. 

[63] M. Mikolas, Farey series and their connection with the prime number problem I. Acta 
Univ. Szeged. Sect. Sci. Math. 13 (1949), 93-117. 

[64] M. Mikolas, Un Theoreme d'equivlence et ses applications, Norske. Vid. Selsk. Forh. 
Trondheim 22 (1950), no. 28, 128-131. 

[65] M. Mikolas, Farey series and their connection with the prime number problem IL Acta 
Univ. Szeged. Sect. Sci. Math. 14 (1951), 5-21. 

[66] L. Mirsky, A theorem on representation of integers in the scale of r, Scripta Math. 15 
(1949), 11-12. 

[67] T. Okada, T. Sekiguchi and Y. Shiota, Applications of binomial measures to power sums 
of digital sums, J. Number Theory 52 (1995), no. 2, 256-266. 

[68] T. Okada, T. Sekiguchi and Y. Shiota, A generalization of Hata-Yamaguti's results on 
the Takagi function II. Multinomial case, Japan J. Indust. Appl. Math. 13 (1996), No. 3, 
435-463. 

[69] Y. Peres, W. Schlag and B. Solomyak, Sixty years of Bernoulli convolutions, in: Frac- 
tal geometry and Stochastics II (Grefswald/Koserov 1998), 39-65, Prog. Prob. Vol. 46, 
Birkhauser: Basel 2000. 

[70] A. Pinkus, Weierstrass and Approximation Theory, J. Approximation Theory 107, 1-66. 

[71] G. de Rham, Sur une courbe plane, J. Math. Pures Appl. Ser. 9 35 (1956), 25-42. 

[72] G. de Rham, Sur quelques courbes definies par des equations fonctionnelles, Univ. e Po- 
litec. Torino. Rend. Sem. Mat. 16 (1956/57), 101-113. 

[73] G. de Rham, Sur un exemple de fonction continue sans deriveee, Enseign. Math. 3 (1957), 
71-72. 

[74] R. Salem, On som,e singular monotonic functions which are strictly increasing. Trans. 
Amer. Math. Soc. 53 (1943), 427-439. 



30 



[75] J. Schauder, Zur Theorie stetiger Abbildungen in Funktionalrdumen, Math. Z. 26 (1927), 
47-65. 

[76] T. Sekiguchi and Y. Shiota, A generalization of Hata-Yamaguti's results on the Takagi 
function, Japan J. Indust. Appl. Math. 8 (1991), 203-219. 

[77] A. Shidfar and K. Sabetfakhri, On the continuity of van der Waerden's function in the 
Holder sense, Amer. Math. Monthly 93 (1986), 375-376. 

[78] A. Shidfar and K. Sabetfakhri, On the Holder continuity of certain functions, Exposition. 
Math. 8 (1990), 365-369. 

[79] Y. Shiota, T. Sekiguchi and T. Okada, Hata-Yamaguti's result on Takagi functions and 
its applications to digital sum problems, In: (Analytic Number Theory , Kyoto 1995) 
(Japanese) , No. 961, Aug. 1996, 73-80. 

[80] K. B. Stolarsky, Power and exponential sums of digital sums related to binomial coefficient 
parity, SIAM J. Appl. Math. 32 (1977), 717-730. 

[81] T. Takagi, A simple example of the continuous function without derivative, from Pro- 
ceedings of the Physico-Mathematical Society of Japan, ser II, Vol 1. 1903, pp 176-177. 
[Collected Papers of Teiji Takagi (S. lyanaga, Ed), Springer Verlag, New York 1990]. 

[82] C. Tricot, Function norms and fractal dimension, SIAM J. Math. Anal. 28 (1997), no. 1, 
189-212. 

[83] J. R. Trollope, An explicit expression for binary digital sums, Math. Mag. 41 (1968), 
21-25. 

[84] L. de Vito, Su un esempio di funzione continua senza derivata, Enseign. Math. 4 (1958), 
53-57. 

[85] B. L. van der Waerden, Ein einf aches Beispiel einer nichtdifferenzierbaren stetigen Funk- 
tion. Math. Z. 32 (1930), 474-475. 

[86] A. Wintncr, Asymptotic distributions and infinite convolutions. Lecture notes by Ralph P. 
Boas, Frank Smithies and John W. Tukey with sympathetic encouragement from Cyrus 
C. Mac Duffcc. The Institute for Advanced Study., Princeton NJ 1938. 

[87] Y. Yamaguchi, K. Tanikawa and N. Mishima, Fractal basin boundary in dynamical sys- 
tems and the Weierstrass- Takagi functions, Physics Letters A 128 (1988), No. 9, 470-478. 

[88] M. Yamaguti and M. Hata, Weierstrass's function and chaos, Hokkaido Math. J. 12 (1983), 
333-342. 

[89] M. Yamaguti, Schauder expansion by some quadratic base function, J. Fac. Sci. Univ. 
Tokyl, Sect lA Math. 36 (1989), no. 2, 187-191. 

[90] M. Yamaguti, M. Hata and J. Kigami, Mathematics of Fractals, Translations of Math- 
ematical Monographs, No. 167, American Mathematical Society: Providence, RI 1997. 
(Translation from Japanese of 1993 original.) 



31 



[91] M. Yamaguti and J. Kigami, Some remarks on Dirichlet problems of Poisson equation, 
Analyse mathematique et applications (J.-L. Lion's 50tli birthday), Gauthier-Villars, 
Montrouge 1988, pp. 465-471. 

Jeffrey C. Lagarias 
Dept. of Mathematics 
The University of Michigan 
Ann Arbor, MI 48109-1043 
email: lagariasSumich . edu 



32 



